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This  paper  addresses  the  design  task  of  characterizing  and  locating  Input/output 
stations  within  a  facilities  layout.  Depending  upon  the  particular  case,  an  Input/output 
station  could  be  a  door  to  an  office  or  a  department,  a  floor  location  where  unit  loads 
enter  and  leave  a  department,  or  an  automated  transfer  station  Integrated  to 
autonomous  robotic  vehicle  systems.  In  all  cases,  the  locations  of  Input/output  stations 
have  a  strong  Impact  on  both  the  cost  of  flow  among  departments  and  the  Internal 
configuration  of  the  departments. 

"  In  the  past,  characterizing  and  locating  Input/output *^1/0)  stations  has  generally 
been  restricted  to  a  secondary  consideration  In  the  detailed  phase  of  facilities  design. 
However,  practical  relevancy  has  recently  lead  researchers  to  more  fundamental 
consideration  of  Input/output  station  locations  when  evaluating  layouts. 

O’Brien  and  Abdel  Barr  (Ref.  9)  have  proposed  a  layout  Improvement  algorithm 
named  S-ZAKY,  similar  to  CRAFT  (Ref.  1),  which  computes  the  expected  flow  distance 
savings  at  each  Iteration  based  on  the  location  of  the  Input/output  (I/O)  stations. 
They  argue  that  this  represents  a  more  realistic  Improvement  assessment  than 
procedures  which  use  Intercentroid  travel  (e.g.  CRAFT).  They  require  an  Input  and  an 
output  at  fixed  relative  locations  for  each  department.  This  has  the  Impact  of 
specifying  much  of  the  Internal  configuration  of  the  department.  Hence,  the  I/O 
station  locations  were  prespeclfled,  used  for  evaluation,  and  not  subject  to  alteration  '  ->r 
optimization  sake.  Warnecke  et  al.  (Ref.  n)  confirmed  that  Interstation  flow  distance 
measurements  are  more  representative  of  actual  flows  than  centroid  to  centroid 
approximations. 

Warnecke  and  Dangelmaler  (Ref.  13)  Incorporated  Into  the  INTALA  Interactive 
layout  software  a  layout  construction  algorithm  permitting  three  Internal  configurations 
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for  each  department.  Each  of  these  three  configurations  fixes  an  input  station  and  an 
output  station  at  given  relative  locations  within  the  department.  At  each  iteration  the 
Interactive  heuristic  selects  a  new  department  to  enter  the  layout,  then  looks  for  the 
best  combination  locatlon/conflguratlon  for  the  department  based  on  the  currently 
located  departments.  Hence,  INTALA  provides  primitive  capability  for  Integrating  I/O 
station  location  and  facility  layout. 


2.  STATION  CHARACTERI2ATION  AND  LOCATION  FOR  A  GIVEN 
LAYOUT 

This  paper  attempts  to  provide  a  systematic  methodology  for  accomplishing  the 
design  task  of  characterizing  and  locating  I/O  stations  for  a  given  facility  layout,  and 
to  demonstrate  how  it  can  potentially  have  a  strong  impact  on  the  overall  flow 
efficiency  of  a  facility.  The  proposed  methodology  Is  summarized  below. 


The  first  task  Is  to  define  the  station  set  and  characterize  each  station.  Simply 
put.  It  means  that  for  each  department  the  designer  must  decide  how  many  stations  are 
to  be  used  In  Its  operation  and  what  flow  types  will  be  assigned  to  each  station.  For 
example,  the  engineer  may  decide  that  a  given  assembly  department  should  work  with 
two  stations,  one  global  Input  station  and  one  global  output  station. 

The  second  task  Is  to  specify  all  interstation  flows  based  on  the  station 
characterizations.  This  Is  similar  to  what  is  commonly  done  In  specifying  department 
flows  for  facilities  layout  studies,  except  that  It  Is  performed  at  the  station  level  Instead 
of  the  department  level.  The  result  Is  a  set  of  Interstations  flows  (e.g.,  500  trlps/perlod 
from  the  output  station  of  the  cell  XYZ  to  the  Input  station  of  storage  zone  ABC). 


As  a  third  task,  the  designer  must  specify  the  location  boundary  region  for  each 
station.  This  Is  usually  based  on  considerations  such  as  technological  feasibility  and 
ease  of  Implementation.  One  extreme  Is  to  Qx  a  station  at  a  given  location  within  the 
department.  At  the  other  extreme,  a  station  may  be  left  free  to  be  located  anywhere 
within  the  department. 

The  fourth  task  Is  to  actually  optimize  the  station  locations  given  the  Hows  and 
boundary  regions.  The  objective  Is  to  obtain  the  set  of  locations  which  maximizes  the 
overall  flow  efficiency  of  the  layout.  .An  efficient  model  and  algorithm  for  this  problem 
is  presented  In  Section  8. 

A  Anal  task  Is  to  perform  sensitivity  analysis  on  the  optimization  results.  The 
purpose  of  the  sensitivity  analysis  Is  to  evaluate  the  robustness  of  the  solution  with 
respect  to  decisions  made  within  the  first  three  tasks. 

3.  STATION  SET  CHARACTERIZATION 

The  first  task  the  designer  has  to  perform  Is  to  determine  the  set  of  stations  to  be 
located,  as  well  as  a  specific  characterization  of  each  station,  which  Is  logical  and 
representative  of  the  system  operation.  In  order  to  facilitate  this  operation.  It  is 
proposed  that  each  department  be  studied  In  order  to  categorize  Its  stations  according 
to  the  following  terminology. 

A  station  will  first  be  categorized  as  either  a  FREE  station  or  a  DEPARTMENT 
station.  A  DEP.ARTMENT  station  belongs  to  a  specific  department  and  only  handles 
flows  associated  with  Its  department.  A  FREE  station  does  not  belong  to  any 
department  and  may  handle  flows  associated  with  any  departments  or  other  stations. 
Instances  where  defining  a  FREE  station  Is  useful  Include  representation  of  Interfloor 
elevators  or  stairs,  Interbulldlng  conveyor  f*ntrances/ex!ts.  critical  doors  or  interfaces 
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within  the  facility,  and  Intersection  nodes  within  the  flow  structure,  a 
DEPARTMENT  station  will  further  be  categorized  according  to  the  types  of  flow  it 
handles.  In  this  development,  a  DEPARTMENT  station  will  be  designated  as  an 
INPUT,  OUTPUT,  INPUT/OUTPUT,  or  INTERNAL  station. 

An  INPUT  station  handles  flow  coming  Into  a  department.  An  OUTPUT  station 
handles  flow  going  out  of  a  department.  An  INPUT/OUTPUT  (I/O)  station  handles 
both  Input  and  output  flows.  Finally  an  INTERNAL  station  has  flows  only  with 
DEPARTMENT  stations  belonging  to  Its  own  department.  Defining  INTERNAL 
stations  Is  especially  useful  with  process  oriented  departments.  For  example,  if  three 
different  sets  of  equipment  are  located  within  a  department  and  there  are  distinct  flows 
associated  with  each  group.  It  Is  frequently  useful  to  define  an  INTERNAL  station 
corresponding  to  each  equipment  set. 

It  Is  also  convenient  to  define  GLOBAL  and  FLOW  GROUP  stations.  A  GLOBAL 
station  is  a  DEPARTMENT  station  through  which  are  directed  all  flow  of  a  given  type 
(l.e..  INPUT.  OUTPUT,  or  INPUT/OUTPUT)  for  Its  department.  For  example,  a 
GLOBAL  OUTPUT  station  Is  a  DEPARTMENT  station  which  deals  with  all  outgoing 
flows  from  Its  department.  A  FLOW  GROUP  station  Is  a  DEPARTMENT  station 
dealing  with  flows  Identified  as  part  of  Its  associated  group.  A  station  which  handles 
only  flows  of  tools  Is  an  example  of  a  FLOW  GROUP  station.  Hence,  GLOBAL  and 
FLOW  GROLT*  stations  are  Identified  by  the  set  of  flows  they  deal  with  within  a  given 
flow  type. 

Finally,  It  Is  convenient  to  designate  stations  In  terms  of  logical  assumptions  made 
about  their  location.  A  station  whose  location  is  restricted  to  be  within  a  boundary 
region  to  be  subsequently  established  will  be  called  a  BOUNDARY  REGION  station. 
In  subsequent  discussions  this  Is  the  default  categorization.  Stations  which  are  assumed 


co  be  at  the  centroid  of  a  department  will  be  called  CENTROID  stations.  FDQSD 
stations  will  refer  to  stations  proposed  to  be  fixed  to  a  specific  location,  different  from 
the  department  centroid. 

In  order  to  Illustrate  the  overall  methodology,  an  example  case  will  be  discussed 
throughout  the  paper.  This  case  Is  adapted  from  Ref.  2.  Relevant  data  Including 
department  descriptions  and  areas  as  well  as  the  Interdepartment  directed  flows  Is  given 
in  Table  l.  Figure  l  presents  the  block  layout  for  the  facility.  Neither  the  building 
shape  nor  the  department  block  layout  have  been  optimized. 

Table  2  summarizes  the  resulting  station  set  characterization.  The  Receiving  and 
Shipping  departments  have  been  attributed  a  single  station,  respectively  GLOBAL 
OUTPUT  and  GLOBAL  INPUT.  The  Milling,  Presses.  Lathes.  Drills,  Welding  and 
Grinding  departments  are  typical  process  departments.  The  multiplicity  of  equipments 
within  each  has  led  the  engineer  to  assign  a  GLOBAL  INPUT  station  and  a  GLOBAL 
OUTPUT  station  to  each  of  these  departments,  coupled  with  a  CENTROID 
INTERNAL  station.  The  flow  logic  Is  as  follows.  Material  arrives  through  the 
GLOBAL  INPUT  station,  then  It  Is  directed  toward  one  of  the  various  equipment  cells. 
Once  the  work  Is  completed  In  a  cell,  the  material  Is  sent  to  the  GLOBAL  OUTPUT 
station  for  final  departure.  Flow  toward  and  from  the  various  cells  Is  approximated  by 
flow  to  and  from  the  CENTROID  INTERNAL  station.  The  Plating  and  Assembly 
departments  can  be  internally  laid  out  to  take  advantage  of  station  locations,  hence 
they  have  been  assigned  GLOBAL  INPUT  and  GLOBAL  OUTPUT  stations.  Since  the 
Internal  configuration  of  the  Stores  and  Warehouse  departments  is  extremely  flexible, 
the  engineer  has  decided  to  permit  FLOW  GROUP  stations  based  on  geographical 
locations  for  these  two  departments.  For  example,  the  Warehouse  department  has  two 
FLOW  GROUP  INPUT  stations.  One  deals  with  Incoming  flows  from  Welding  and 
Plating  while  the  other  deals  with  flow  from  .Assembly. 
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TABLE  2 


STATION  CHARACTERIZATION  AND  BOUNDARY  SPECIFICATION 


STATION 


CHARACTERISTICS 


Department 

# 

Flow 

Flow 

Location 

K 

Type 

Set 

Assumption 

fe 

Receiving 

1 

O 

GL 

B 

Milling 

1 

I 

GL 

B 

2 

INT 

GL 

C 

5 

3 

O 

GL 

B 

s 

Presses 

1 

I 

GL 

B 

2 

INT 

GL 

C 

3 

O 

GL 

B 

Us 

Lathe* 

1 

I 

GL 

B 

3 

INT 

GL 

C 

3 

O 

GL 

B 

Drills 

1 

I 

GL 

B 

£ 

2 

INT 

GL 

C 

£ 

3 

O 

GL 

B 

Welding 

1 

I 

GL 

B 

2 

ENT 

GL 

C 

3 

O 

GL 

B 

e 

Plating 

1 

I 

GL 

B 

2 

O 

GL 

B 

Grinding 

1 

I 

GL 

B 

2 

INT 

GL 

C 

£ 

3 

O 

GL 

B 

Assembly 

1 

I 

GL 

B 

3 

O 

GL 

B 

i 

Warehouse 

l 

I 

FG 

B 

2 

I 

FG 

B 

Shipping 

1 

I 

GL 

B 

K 

Stores 

1 

I 

GL 

B 

3 

O 

FG 

B 

'•K 

3 

O 

FG 

B 

BOUNDARIES 
X-Wil  Y-» 

Min.  Mix 


0.7* 

0,53 

10.6.18.5 
0.33 
33,30 
01.01 
33,00 
0.100 
SO, 60 
0,100 

160.300 

183.6.103.6 

160.300 
80.110 

103.6.103.6 
80.110 
100.160 
110.160 
110.174 
146.146 
110.174 
110.174 
110.174 
78.110 


117.140 
30.01 

64.64 
30.01 
30.73 

64.64 
30.73 
0,30 
10.18 
0.30 
0.01 

30.6.30.6 
0.01 
30.73 

64.64 
30,72 
0.33 
0.72 

106.140 

123.6.133.6 

106.140 
72.106 
72.106 

72.140 


FLOW  GROUP 
Specification 


Prom  Wcldlnt 
ft  Pitting 
From  Assembly 


To  Press 
To  Lathes 


ABBREVIATIONS 

FLOW  TYPE:  I— INPUT.  O— OUTPUT.  INT— INTERNAL 
FLOW  SET:  GL— GLOBAL,  FG— FLOW  GROUP 
LOCATION  ASSUMPTION:  B— BOUNDARY.  C- CENTROID 


4.  FLOW  SPECIFICATION 


Once  all  stations  to  be  located  have  been  characterized,  the  designer  must  specify 
the  magnitude  of  all  Interstation  flows.  This  process  Is  similar  to  specification  of  the 
standard  Interdepartment  FROM-TO  and  FLOW-BETWEEN  charts  (for  background 
reference,  see  Ref.  12).  For  example,  rather  than  determining  flow  from  a  department 
A  to  a  department  B  for  a  FROM-TO  chart,  flow  Is  determined  from  each  OUTPUT 
station  of  department  A  to  each  INPUT  station  of  department  B.  If  the  designer  limits 
himself  to  a  single  GLOBAL  I/O  station  for  every  department,  then  the  flow  list  will 
correspond  directly  to  the  standard  lnterdepartment  FLOW-BETWEEN  chart. 
Furthermore.  If  he  defines  two  stations  for  every  department,  one  as  a  GLOBAL 
INPUT  station  and  one  as  a  GLOBAL  OUTPUT  station,  then  the  flow  list  win 
correspond  directly  to  the  standard  lnterdepartment  FROM-TO  chart. 

Table  3  presents  the  resulting  flow  set  for  the  example,  based  on  the  station 
characteristics  expressed  In  Table  2  and  the  lnterdepartment  FROM-TO  flows  from 
Table  1.  The  lnterstatlon  flows  are  flows  between  stations  of  distinct  departments. 
These  flows  are  straightforward  to  determine.  For  example.  Table  1  shows  that  there  Is 
a  flow  of  30  from  Receiving  to  Milling.  Receiving  has  a  single  GLOBAL  OUTPUT 
station  (RE,  1),  and  Milling  has  a  GLOBAL  INPUT  station  (MI,  1),  hence  there  Is  a 
flow  of  30  from  station  (RE,  1)  to  station  (MI.l  ).  As  another  example,  there  Is  a  flow 
of  50  from  .Assembly  to  Warehouse.  Warehouse  has  a  FLOW  GROUP  INPUT  station 
(WR,  2)  specialized  In  dealing  with  flow  Incoming  from  .Assembly.  Assembly  has  a 
GLOBAL  OUTPUT  station  (AS,  2).  Therefore,  there  Is  a  flow  of  50  from  (AS.  2)  to 
(WR.  2). 

The  lntradepartmental  flows  are  flows  between  stations  within  a  given  department. 
Table  2  specifies  that  the  Lathes  department  has  GLOBAL  INPUT,  INTERNAL  and 
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TABLE  3 

interstation  flow  set  specification 


1 

INTERSTATION  FLOWS 

INTERDEPARTMENT 

INTRADEPARTMENT 

STATION 

station 

FLOW 

STATION 

STATION 

FLOW 

RE.l 

MU 

30 

MI.l 

M1.3 

30 

Isi 

RE.l 

ST.l 

30 

MI.2 

MLS 

80 

MI. 3 

PR.l 

40 

PR.l 

PR.3 

80 

MI. 3 

LA.l 

10 

PR.3 

PR.S 

80 

g 

MI. 3 

PL.l 

10 

LA.l 

LA.3 

S6 

PR.3 

LA.l 

36 

LA. 2 

LA-3 

76 

PR.3 

DR.l 

6 

DR.l 

DR.3 

36 

& 

LA. 3 

DR.l 

20 

DR.3 

DR.3 

36 

LA  .3 

WE.l 

10 

WE.l 

WE.3 

36 

LA, 3 

PL.l 

46 

WE.3 

WE.3 

36 

& 

DR.3 

PL.l 

S 

GR.1 

GR.2 

86 

■£ 

DR.3 

GR.1 

20 

GR.2 

GR.3 

106 

Dr. 3 

AS.l 

10 

y 

WE  .3 

LA.l 

10 

£ 

>  WE.3 

AS.l 

16 

WE. 3 

WR.l 

10 

PL.2 

WE.l 

6 

*.* 

PL.2 

GR.1 

46 

. j 

PL.2 

AS.l 

10 

PL.2 

WR.l 

20 

ft 

GR.3 

WE.l 

20 

1 

1  GR.3 

AS.l 

26 

GR.3 

SH.l 

80 

V*, 

AS.2 

WR.2 

60 

V, 

20 

V. 

ST.2 

PR.l 

►  k> 

ST.3 

LA.l 

30 

p 


; 


v 

%■ 


OUTPUT  stations.  Given  the  flow  logic  expressed  earlier,  all  flow  coming  Into  the 
Lathes  department  through  (LA,  1)  Is  approximated  by  flow  to  the  INTERNAL  station 
(LA,  2).  By  summing  all  flows  to  (LA,  l),  the  flow  from  (LA,  1)  to  (LA,  2)  can  be 
computed  as  85.  The  same  logic  leads  to  a  flow  of  75  from  (LA,  2)  to  the  GLOBAL 
OUTPUT  station  (LA,  3).  Material  transformation,  scrap,  unit  load,  lot  sizing, 
handling  methods  may  all  Justify  the  difference  between  the  total  Input  and  output 
flows. 

5.  BOUNDARY  REGION  SPECIFICATION 

It  Is  generally  desirable  to  specify  a  boundary  region  within  which  a  station  can  be 
located.  .An  infinite  boundary  region  for  a  station  means  that  the  station  can  be 
located  anywhere  within  the  facility.  A  point  boundary  region  indicates  that  the 
station  Is  to  be  fixed  at  a  prespeclfled  location.  For  example,  a  CENTROID  station  is 
assumed  to  be  fixed  at  a  point  which  Is  the  department  centroid.  Departments  where 
the  input  and  output  are  technologically  predetermined,  such  as  product  line  assembly 
departments  when  the  layout  Is  difficult  to  change,  may  also  have  a  point  boundary 
region. 

Department  boundaries  are  typically  used  as  station  boundaries  when  the  designer 
does  not  conceive  major  restrictions  on  station  locations.  Also  In  very  aggregate  studies, 
the  station  boundaries  can  usually  be  taken  as  the  department  boundaries.  In  less 
aggregate  studies,  the  designer  may  compensate  for  technological  constraints  on 
stations  by  restricting  them  to  be  a  given  distance  within  the  department  boundary 
region. 

For  convenience,  the  resulting  station  boundaries  for  the  example  case  have  been 
Included  In  Table  2.  All  boundary  regions  are  presented  as  the  extreme  X-axls  and  Y- 


axis  coordinates  of  a  rectangle.  The  limitation  to  rectangular  boundary  regions  will  be 
Justified  In  a  latter  section.  The  CENTROID  INTERNAL  stations  have  point 
boundary  regions  at  their  department  centroid.  The  designer,  not  conceiving  a  priori 
any  major  restriction  on  location  for  other  stations,  has  used  the  department 
boundaries  as  boundary  regions.  This  can  be  visualized  by  referring  to  the  layout  in 
Figure  1. 

0.  STATION  LOCATION  OPTIMIZATION 

Given  a  set  of  stations,  each  restricted  to  be  within  a  specified  boundary  region, 
the  location  optimization  task  Is  to  find  the  set  of  station  locations  that  permit  the 
most  effective  Interstation  flows. 

The  best  measure  of  efficiency  to  be  used  for  station  location  Is  subject  to  debate. 
Extremely  detailed  measures  of  efficiency  can  be  obtained  by  simulation  studies  of  a 
given  set  of  station  locations.  However,  In  order  to  determine  the  locations.  It  Is 
necessary  to  use  a  more  aggregate  measure  of  efficiency. 

The  optimization  criterion  proposed  here  Is  to  minimize  the  sum  of  rectilinear  (Ll) 
distance  travelled  by  all  interstation  flows.  It  provides  a  reasonable  approximation  to 
the  actual  objective  and  still  yields  a  tractable  model.  Under  the  assumptions  of  a 
convex  piece-wise  linear  boundary  region  for  each  station,  the  location  problem  can  be 
solved  using  a  linear  programming  model.  Under  the  additional  assumption  of  a 
rectangular  boundary  region  for  each  station,  a  much  more  efficient  network  theory 
based  solution  methodology  can  be  applied.  The  Impact  of  these  assumptions  on  model 
realism  will  be  discussed  subsequently. 

The  linear  programming  mode!  for  optimizing  the  station  locations  under  Ll 
criterion  and  with  rectangular  boundary  regions  Is  presented  below.  First  let 
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The  station  location  problem  can  then  be  stated  as: 
Minimize 
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2k  <  xi  <  x, 

v{ie  s} 

(3) 

Xi  <  Yi  <  Vi 

v{i  e  s} 

Minimizing  expression  1  corresponds  to  minimizing  the  sum,  over  all  stations  having  a 
positive  flow,  of  the  product  of  the  rectilinear  (Ll)  distance  between  the  stations  and 
their  flow  value.  It  will  be  referred  to  as  the  Ll  flow  score.  Expressions  2  and  3 
guarantee  that  the  stations  are  located  within  their  boundary  rectangular  grid  region. 
It  Is  a  simple  matter  (Ref.  3)  to  reformulate  expression  1  as  a  linear  function  with 
additional  variables  and  linear  constraints. 

Expressions  2  and  3  can  be  altered  to  obtain  a  linear  programming  formulation  as 
long  as  the  station  locations  are  restricted  to  be  within  any  specified  convex 
polyhedron.  As  a  simple  numerical  example,  station  1  can  be  restricted  to  be  located 
within  the  triangle  formed  by  points  (5,8)  (10,12)  and  (15,4)  by  using  the  following 
expressions  defining  the  three  line  equations  forming  the  triangle  Internal  region: 

y,  -  (8/5)  x,  <  0 

y,  +  (1/5)  x,  >  7 
y,  +  (8/5)  x,  <  28 
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7.  RELATIONSHIP  WITH  MULTIFACILITY  LOCATION  MODEL 


The  model  expressed  by  (1)  to  (3)  Is  equivalent  to  the  well-known  and  very 
efficiently  solved  rectilinear  multlfaclllty  location  problem  (Ref.  11).  The  rectilinear 
multlfaclllty  location  problem  Is  defined  as  follows.  A  set  of  new  facilities  ,  with  no 
location  constraints,  are  to  be  located  In  order  to  minimize  the  sum  of  the  rectilinear 
distance  travelled  by  all  flows  between  the  new  facilities,  and  between  the  new  facilities 
and  a  set  of  existing  facilities  specified  by  given  locations. 

It  Is  easy  to  see  that  the  rectilinear  multlfaclllty  model  is  a  special  case  of  the 
station  location  model.  The  set  of  new  facilities  corresponds  to  a  set  of  FREE  stations 
with  infinite  boundary  regions.  The  set  of  existing  facilities  corresponds  to  a  set  of 
FREE  stations  with  point  boundary  regions. 

To  see  that  the  models  are  In  fact  equivalent  let  any  station  with  a  point  boundary 
region  correspond  to  an  existing  facility,  and  any  station  with  an  infinite  boundary 
region  correspond  to  a  new  facility.  The  complication  occurs  when  stations  have  a 
rectangular  boundary  region.  To  enforce  the  boundary  region  an  existing  facility  can 
be  created  at  each  of  the  four  corners  of  the  rectangular  boundary  region  with  a  large 
(BIG-M)  flow  defined  between  the  new  facility  and  each  of  the  four  created  existing 
facilities.  Using  the  properties  derived  by  Picard  and  Ratliff  (Ref.  11),  locating  the  new 
facility  anywhere  within  the  rectangular  region  between  the  four  existing  facilities  leads 
to  equivalent  optimal  solutions  for  the  BIG-M  flows  (score  related  to  these  Is  identical). 
Furthermore,  the  new  facility  cannot  be  located  outside  the  4-faclIlty  region  in  a  global 
optimal  solution  since  the  BIG-M  flows  dominate  the  original  flows.  Hence  the  stations 
are  restricted  to  their  rectangular  boundaries  within  the  framework  required  by  the 
rectilinear  multlfaclllty  location  model. 


S.  OPTIMIZATION  METHODOLOGY  WITH  RECTANGULAR  BOUNDARIES 
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The  transformation  process  described  In  section  7  allows  development  of  a 
specialized  algorithm  for  solving  station  location  optimization  under  Ll  criterion  with 
rectangular  boundaries.  The  algorithm  Is  a  variation  of  the  algorithm  designed  by 
Picard  and  Ratliff  (Ref.  11)  for  rectilinear  multlfaclllty  location.  The  proposed 
algorithm  for  station  location  Is  presented  In  figure  2. 

The  location  optimizations  for  the  horizontal  and  vertical  dimensions  are 
performed  Independently.  For  each  dimension,  the  algorithm  moves  Iteratively  from 
the  lowest  boundary  coordinate  to  the  highest  boundary  coordinate.  At  each  boundary 
coordinate  C,  It  divides  the  stations  among  three  sets.  The  set  «—  includes  all  stations 
whose  upper  boundary  Is  smaller  than  C,  these  are  already  located  due  to  the  increasing 
order  of  treatment  of  C.  The  set  — ►  Includes  all  stations  whose  lower  bound  Is  larger 
than  C,  these  cannot  be  located  at  C.  The  set  ©  Includes  all  other  stations,  in  fact  all 
those  which  may  be  located  at  C. 

Given  these  three  station  sets,  the  algorithm  creates  a  graph  G  as  follows.  A  node 
1  Is  defined  to  Include  all  stations  within  set  «—  and  a  node  2  Is  defined  to  Include  all 
stations  within  set  — Then  a  node  Is  defined  to  represent  each  station  within  set  ©. 
A  link  Is  then  defined  between  each  pair  of  nodes  (m  and  n).  The  capacity  of  a  link, 
Cmn  ,  Is  set  equal  to  the  sum  of  the  flows  (f,,)  between  all  stations  Included  in  node  m 
and  all  stations  Included  In  node  n. 

The  algorithm  then  finds  the  MINIMUM  CUT  (Ref.  10)  between  nodes  l  and  2  In 
graph  G,  based  on  the  link  capacities  (Cmn)  .  The  cut  Is  denoted  as  L/L  ,  with  node 
1  €  L  .  The  cut  divides  the  nodes  In  two  sets  L  and  L.  Hence  the  stations  are  divided 
In  two  sets  L  and  L  .  The  set  L  Includes  all  stations  In  set  *—  and  stations  In  set  © 
which  have  more  flow  with  the  other  stations  In  L  than  with  all  those  In  L.  A  station 
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BEGIN 

FOR  EACH  DIMENSION  INDEPENDENTLY  (hor.  and  vert.)  DO 
BEGIN 
REPEAT 

FOR  EACH  BOUNDARY  COORDINATE  C.  IN  INCREASING  ORDER.  DO 
BEGIN 

1.  LOCATE  AT  C  ALL  UNLOCATED  STATIONS  WHOSE  UPPER  BOUNDARY 

2.  DIVIDE  ALL  S  STATIONS  AMONG  3  SETS  —  AS  FOLLOWS: 

*— :  {stations  whose  upper  boundary  <  C} 

©:  {stations  whose  lower  boundary  <  C 
and  whose  upper  boundary  >  C} 

— {stations  whose  lower  boundary  >  C} 

3.  DEFINE  A  GRAPH  G  AS  FOLLOWS: 

•  node  1 :  {Stations  €  *— } 

•  node  2:  {Stations  6  — •  } 

•  node  (3 . [|©  ||  t 2)  for  each  station  £© 

•  a  link  between  each  pair  of  nodes  (m  and  Q), 
with  capacity  (cmn)  equal  to  the  sum  of  the 
flows  (f(J)  between  all  stations  16  m 

and  J  €  n. 

4.  FIND  THE  MINIMUM  CUT  (REF.  10)  BETWEEN  NODES  1  AND  2  IN 
GRAPH  G.  BASED  ON  THE  LINK  CAPACITIES  (cmn): 

DENOTE  THE  CUT  BY  L/L  .  WITH  NODE  I  6  L 

5.  LOCATE  AT  C  ALL  STATIONS  1  €©  WHOSE  NODE  (1+2)  £L 
AND  FDC  THEIR  BOUNDARIES  AT  C 

END 

UNTIL  ALL  S  STATIONS  ARE  LOCATED 
END 
END 


FIGURE  2  Proposed  alcorithm  for  the  rectilinear  station  location 

model  with  end  rectangular  boundary  regions 


!n  set  ©  which  ends  up  In  L  Is  more  attracted  to  the  lower  coordinate  locations  than 
hlgher-than-C  coordinate  locations.  However,  this  was  not  true  at  the  previous  C 
Investigated,  otherwise  It  would  already  be  located.  Hence  this  Indicates  that  such  a 
station  should  be  located  exactly  at  C.  For  algorithmic  purposes  Its  boundaries  are 
then  fixed  at  C.  The  algorithm  then  moves  to  the  following  C  and  repeats  the  process 
until  all  stations  are  located. 

Figure  3  presents  an  application  of  the  methodology  for  the  example  case  under 
the  simpler  assumption  of  a  single  GLOBAL  I/O  station  for  each  department.  The 
lnterstatlon  flows  f„  are  non-dlrectlonal  (between),  so  flow  from  1  to  J  and  from  J  to  1 
are  added  to  get  f,,  .  The  boundary  region  defined  for  each  department's  station  Is 
depicted  graphically  by  an  horizontal  boundary  line  and  a  vertical  boundary  line.  Both 
horizontal  and  vertical  optimizations  are  detailed.  For  each  Iteration,  the  three  sets  , 
©  and  —  are  listed,  then  the  minimum  cut  computed  Is  Included,  Anally  the  resulting 
location  decisions  are  shown.  The  final  optimal  locations  are  graphically  represented  on 
the  layout  and  the  associated  Ll  flow  score  Is  tabulated. 

Cases  with  over  100  stations  have  been  solved  within  100  seconds  on  a  IBM  AT 
personal  computer  with  a  mathematical  co-processor,  and  within  a  few  seconds  on  a 
CDC  CYBER  17-4.  Hence,  the  methodology  can  be  applied  to  most  Industrial  design 
cases  without  bothering  about  optimization  solution  time  within  the  overall  design 
process. 

9.  SCORE  BOUNDING 

As  part  of  the  suggested  sensitivity  analysis  phase.  It  may  be  worthwhile  to 
compute  lower  and  upper  bounds  on  the  flow  score  given  the  flows  and  boundary 
regions  In  order  to  assess  how  restrictive  are  the  Imposed  constraints. 
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A  lower  bound  on  the  Ll  flow  score  can  be  obtained  by  assuming  for  each  pair  of 
stations  having  a  positive  flow  that  they  are  located  so  that  the  distance  between  them 
is  as  small  as  possible.  Summing  the  obtained  distance  times  the  flow  over  ail  pairs 
Independently  gives  a  lower  bound  Su  on  tbe  Ll  score.  Taking  advantage  of  the 
rectangular  boundaries,  Sn  can  be  computed  efficiently  as  follows: 

(4)  £u  =  S  ru  f max  (to  “  *j)  •  ~  *i)  •  °)  +  max(to  "  ?i)  *  ~  ?j)  •  °)j 

F  1  J 

By  setting  the  (&*  ,  xk)  and  (x*  ,  ?k)  equal  to  their  largest  possible  region,  which  usually 
corresponds  to  the  department  boundaries,  £li  permits  to  extend  the  boundary  to  the 
lowest  Ll  score  attainable  with  full  degree  of  freedom  In  station  design  given  the  actual 
layout. 

Inversely,  by  assuming  for  each  pair  of  stations  having  a  positive  flow  that  they 
are  located  as  far  from  each  other  as  possible,  and  then  summing  those  distances  times 
the  flow  over  all  pairs  Independently  gives  an  upper  bound  SLl.  Taking  again 
advantage  of  the  rectangular  boundaries,  can  be  computed  as  follows: 

(5)  Su  =  £  f,»  [  max((x,  -  *,) .  (I,  -  &))  +  max  ((y,  -  y,) ,  (V,  -  y,))  1 

F  L  J 

SLl  corresponds  to  the  worst  possible  station  design  strategy  by  locating  stations  always 
as  far  as  possible  from  each  other.  Determining  this  upper  bound  Is  useful  to  And  out  if 
the  Ll  flow  score  Is  sensitive  upon  the  set  of  station  locations.  If  Su  and  are  not 
significantly  different,  then  optimizing  the  station  location  will  not  have  a  significant 
Impact  on  the  flow  score.  Then  the  actual  station  locations  may  be  based  on  other 
considerations  than  flow. 

For  the  example  of  figure  3,  £Ll  Is  11,950.  and  $Ll  Is  105.505.  The  Ll  flow  score 
obtained  for  the  optimal  solution  Is  20.600.  which  Is  25%  of  the  upper  bound  and  223% 
of  the  lower  bound.  The  score  being  therefore  more  than  twice  the  lower  bound.  It 
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Indicates  that  the  single-station  operating  constraint  Is  strongly  affecting  the  expected 
distance  travelled. 

10.  RECTANGULAR  BOUNDARY  REGION  ASSUMPTION 

On  first  observation,  restricting  the  station  location  boundary  region  to  be 
rectangles  may  appear  very  limiting.  However,  when  analyzing  actual  layouts,  one 
rarely  observes  elaborate  department  shapes.  In  fact,  it  Is  good  layout  practice  to  keep 
the  department  shapes  as  simple  as  possible.  The  majority  of  non-rectangular 
departments  are  composites  of  at  most  three  rectangular  areas.  For  example, 
department  AS  In  Figure  3  can  be  viewed  as  a  composite  at  two  rectangular  areas. 

There  are  several  approaches  to  modeling  the  station  location  constraint  for  these 
nonrectangular  departments  as  rectangles.  Sometimes  It  Is  straightforward  to  define  a 
rectangular  region  Inside  the  department  within  which  the  optimal  location  Is  to  be 
located  with  very  high  probability. 

If  it  Is  not  obvious  how  to  define  the  boundary  region  within  an  interior  rectangle, 
an  enclosing  rectangle  can  be  defined  and  used  In  the  optimization.  If  the  optimum 
location  turns  out  to  be  Inside  the  department  then  the  problem  Is  solved.  If  not,  then 
It  Is  often  possible  to  use  knowledge  of  the  optimum  to  define  an  Interior  rectangle  or 
at  most  a  couple  of  possibilities  which  need  to  be  enumerated. 

It  therefore  appears  that  a  rectangular  boundary  assumption  can  be  used  In 
practice  without  significant  loss  of  applicability. 


£ 


* 


* 

a 


i 


« 

* 

£ 


5 


» 


■wiru  WMir.MT.ji 


11.  RESULTING  DESIGN  FOR  THE  EXAMPLE  CASE 

Given  the  station  characterization  and  boundary  specification  expressed  In  Table  2, 
and  the  Interstation  flow  set  specification  as  stated  In  Table  3,  the  optimization  model 
results  In  the  design  pictured  In  figure  *1  after  a  computing  time  of  52  seconds  on  an 
IBM  PC  AT  with  mathematical  co-processor. 

The  resulting  flow  score  Is  35730,  this  Is  higher  than  the  score  of  26600  obtained 
for  the  simple  single  station  assumption.  This  Increase  Is  directly  related  to  the  use  of 
CENTROID  INTERNAL  stations  to  represent  adequately  Internal  movement  through 
the  pertinent  process  departments,  as  decided  by  the  designer. 

It  Is  Interesting  to  note  that  the  algorithm  has  suggested  that,  for  4  out  of  6  of  the 
process  departments  with  a  CENTROID  INTERNAL  station,  their  stations  should  all 
be  clustered  as  a  centralized  Input/output  location  within  the  department.  This  was 
not  suggested,  however,  for  the  Milling  and  Lathes  departments.  The  design  provides 
Illustration  of  typical  clusterings  of  stations  from  different  departments  along  their 
common  boundary,  hence  providing  minimal  travel  for  their  associated  flows.  An 
example  Is  the  clustering  of  the  OUTPUT  station  of  Lathes  and  the  INPUT  station  of 
Plating. 

One  can  note  In  Table  2  that  the  boundary  regions  for  both  stations  of  Plating 
had  to  be  studied  carefully,  even  If  a  full  freedom  of  location  within  the  department 
was  suggested.  This  Is  due  to  the  L  shape  of  the  department.  The  stated  boundaries 
were  decided  upon  by  studying  the  boundary  region  for  the  stations  having  flow  with 
the  stations  of  the  Plating  department.  As  a  sensitivity  test,  the  engineer  solved  the 
case  with  (PL,1)  restricted  to  X  :  (l  16,1561,  Y  :  [0,72].  The  suggested  location  for 
(PL.l)  was  (116,18)  which  confirms  Its  attraction  to  a  lower  horizontal  coordinate;  as 
expected,  the  flow  score  Increased  to  38530. 
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Finally,  the  algorithm  has  taken  advantage  of  the  added  freedom  allowed  by  the 
designer  when  he  defined  distinct  FLOW  GROUP  stations  for  the  Stores  and 
Warehouse  departments.  These  FLOW  GROUP  stations  ended  up  In  distinct  locations 
In  the  suggested  design. 

12.  CONCLUSION 

This  paper  has  proposed  a  methodology  for  characterizing  and  locating 
Input/output  stations  within  a  facility  layout.  The  methodology  Is  to  first  define  the 
station  set  and  characterize  each  station:  second,  specify  all  Interstation  flows  based  on 
the  station  characterizations:  third,  specify  the  location  boundary  region  for  each 
station:  fourth,  optimize  the  station  locations  given  the  flows  and  boundary  regions; 
and,  finally,  perform  sensitivity  analysis  on  the  optimization  results.  A  specialized 
algorithm  was  proposed  to  solve  the  optimization  problem  of  minimizing  the  rectilinear 
flow  score  given  rectangular  boundary  regions. 

The  material  presented  In  this  paper  has  been  used  through  various  empirical 
research  studies  with  Industrial  firms  (Ref.  5,  0,7,8)  and  has  led  In  all  cases  to 
significant  and  representative  designs.  However,  there  are  a  number  of  areas  which 
seem  productive  for  further  research. 

At  the  optimization  model  level,  further  research  Is  necessary  In  addressing  other 
optimization  criteria,  more  complex  boundary  regions,  and  the  cost  of  relocating 
stations  In  an  actual  layout.  Furthermore,  certain  station  set  definition  and  station 
characterization  decisions  could  be  subject  to  optimization  modeling. 

The  location  methodology  presented  here  minimizes  the  expected  distance 
travelled.  Impacts  of  this  are  resulting  station  clusterings  and  high  flow  links  between 
stations.  These  have  some  potentially  very  Interesting  Impacts  on  factory  handling 
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automation  which  could  be  studied  In  more  depth.  As  a  negative  impact  In  some 
settings.  It  may  create  congestion.  Simulation  based  approaches  to  deal  with  this 
phenomenon  could  be  designed.  Empirical  studies  could  also  be  performed  to  assess 
typical  station  characterizations  in  Industry,  Impacts  of  design  constraints  to  overall 
efficiency  achieved,  and  so  on.  These  would  require  extensive  Industry  collaboration  to 
insure  validity  and  representatlvlty  of  the  results. 

Finally,  as  suggested  In  Ref.  4  and  13,  a  promising  avenue  for  further  research  Is 
the  development  of  methodologies  and  models  to  sustain  the  Integrated  design  of 
facility  layout  and  station  location.  The  methodology  presented  In  this  paper  could  be 
used  as  a  building  block  for  this  research. 
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